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1. Introduction 

Log and loglog corrections to physical quantities in critical phenomena generally appear 
in a statistical system at the critical dimension. In the language of renormalization 
group (RG) [1], those corrections arise from marginally irrelevant coupling constants 
in the system. As the length scale we are looking at becomes larger, the coupling 
constants effectively change, obeying a renormalization-group equation (RGE), and 
approach an infrared fixed point if initial values of the trajectories are on the critical 
surface. Universality (i.e., property independent of the initial values) of the logarithmic 
corrections is closely related to the long-distance behavior of the running coupling 
constants. 

Generally, the log and loglog corrections are obtained respectively from the leading 
and the next-to-leading order of the beta function specifying the RGE. In order to 
see this, it is instructive to consider the case of a single coupling constant. When the 
origin is an infrared fixed point, an RGE of a marginal coupling constant g up to the 
next-to-leading order is generally described as 

^L = -ag 2 + bg 3 , a > 0, (1) 

where t is related to a length scale L of RG transformation (RGT) by t — InL. Here 
we consider the weak-coupling region between the two fixed points, < g < a/b. It is 
readily integrated and the solution g(t) satisfies 



at + C = — — In 

g(t) a 



g(t) (£-*(*))" 



(2) 



in this region. The constant C is determined from an initial condition. When t — > oo, 
g(t) has the following asymptotic form: 

/ x 1 bint 

which implies that the first term contributes to a InL correction, while the second 
term to a In InL correction. A key feature is that the coefficients of \jt and \nt/t 2 are 
independent of C, which leads to universal logarithmic corrections. 

Although we can integrate RGE explicitly in the case of a single coupling constant, 
we cannot generally perform the same procedure in the case of multiple coupling 
constants. Therefore, it is worthwhile to determine an asymptotic form analogous to 
(3) in the case of multiple marginal coupling constants, which is the main subject of 
this paper. 

One cannot linearize RGE about a fixed point in the case of marginal coupling 
constants, which complicates the problem finding an asymptotic form without explicit 
integration. An algebraic method was found in Refs [2, 3], where the beta function is 
restricted to the lowest order. Since the lowest-order beta function for marginal coupling 
constants is homogeneous, the RGE is invariant under scaling transformation [2]. One 
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can define another RG transformation to the RGE, thanks to the scale invariance. f 
The new RGE generally has a linear term, which allows us to obtain the asymptotic 
form without explicit integration. 

However, we cannot apply the above method when higher orders of the beta function 
are taken into account because there are no such scale invariances. Hence we need to 
find an alternative method to remedy the problem for linear terms to vanish. 

In the next section, we present a change of variables in the RGE that allows us to 
apply the linearization. In section 3, we switch from the resultant RGE to an equivalent 
integral equation, and outline the existence of a unique solution driven to the fixed 
point. We also give an estimation to the solution. Details of the proof are found in 
Appendix A. In section 4, using the estimation found in the preceding section, we show 
a sufficient condition for loglog corrections be universal. A universal asymptotic formula 
for the solution in the long-distance limit is also derived under the sufficient condition. 
In section 5, applying our result, we rederive the universal asymptotic formula for the 
running coupling constants in the classical XY model, as an example. The result is 
consistent with the original article by Amit et al. [6]. The final section is devoted to 
summary and discussion. 



2. Changing variables of RGE 

We consider an RGE for marginal coupling constants denoted by g(t) = (gi(t), ...,g n (t)). 
We regard the space of the coupling constants as the n dimensional Euclidean space R n . 
Suppose that we have obtained the RGE up to the next-to-leading order, which is to 
say we start with the following RGE 

^ = V(jg(t)) + F(jg(t)). (4) 

The leading and the subleading terms of the beta function are described by V and F 
respectively. It is assumed that they possess the following scaling property: 

V(kg) = k 2 V(g), F(kg) = k 3 F(g). (5) 

In general, V and F are obtained as quadratic and cubic polynomials in the coupling 
constants respectively. The beta function is defined on the whole space R n in this case. 
However, in some cases, the beta function is obtained in a rational form (e.g., see [7]). 
For this reason, it is suitable to assume that V and F are defined on some region E in 
iT.§ 

It is a general feature of an RGE of marginal coupling constants that there are no 
linear terms, which causes difficulty in deriving an asymptotic formula. We introduce 
new variables to bypass this problem. First we replace t by 

u = i log {et + 1) , (6) 

% A general idea of RG, applied as a tool for asymptotic analysis of non-linear differential equations, 
is developed in Refs. [4, 5]. 

§ More precisely, E is open subset of R n whose closure contains the origin. We also assume that V 
and F belong to C 2 (E), i.e., their second derivatives exist and are continuous on E. 
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where e is a parameter with 

< e < 1. (7) 

As we will see later, e is introduced to control an effect of the subleading term F. Next 
we change g to c, where 

g{t) = e-™c(u). (8) 

The left-hand side of (4) becomes 

at du v 

while the right-hand side is 

V(e- £U c) + F(e- tu c) = e~ 2eu V(c) + e~ 3tu F(c) (10) 
because of the scaling property (5). In this way, the RGE (4) is written as|| 

= ec(u) + V(c(u)) + e-™F(c(u)). (11) 

du 

Now we extract the linear part from the first two terms. We assume that a non- 
trivial solution c* G E for 

ec* + V (c*) = (12) 
exists. Note that c* is linear in e. In fact, let us introduce a* by 

c* = ea*. (13) 
Because of the scaling property (5), a* is determined by 

a* + V(a*) =0. (14) 

Hence a* is independent of e and (13) indicates that c* approaches the origin as e 
becomes smaller. Therefore, an effect of the subleading term F in a neighborhood of c* 
is suppressed if we take e sufficiently small. It plays an important role in showing the 
existence of a solution for (18) by a contraction map, as we will see in Appendix A. 
We analyze (11) in a neighborhood of c*. Define 

b(u) = c(u) - c* (15) 

and write 

V{c{u)) = V(c*) + DV{c*)b{u) + v(b(u)), (16) 
where DV(c*) is the derivative of V at c*, which is represented by the n x n matrix as 

W(c%- = g(c*). (17) 

| We comment that this transformation is generalized in the case where a beta function starts with 
an mth homogeneous function (to > 2). In fact, if we define u and c(u) by (m — l)eu = ln(et + 1) and 
g{t) = e~ eu c(u) respectively, we get a similar equation having a linear term. 
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The RGE (11) is written as 
db(u) 



du 

where 



Mb(u) + H(u,b(u)), (18) 



M = el n + DV (c*) , H(u,b(u)) = v(b{u)) + e- tu F(c* + b(u)) (19) 
with J n being the n x n unit matrix. Further, (5) yields^ 

DV(c*) = eDV{a*). (20) 
Therefore M is linear in e 

M = e(/„ + W(a*)). (21) 
To sum up, if we know a solution b(u) for (18), we get a solution for (4) by 

g (t) = e~ eu (c* + 6 (u)) , M=iln(et + 1). (22) 
3. Integral equation 

In this section, following reference [8] (Section 2.7), we derive an integral equation 
satisfied by a solution b(u) for (18) driven to 0. 

We assume that there are no eigenvalues with zero real part in M. Suppose that 
M has k eigenvalues having a negative real part, and n — k eigenvalues with positive 
real part. If the positive eigenmodes are fine-tuned to vanish, \b(u)\ becomes smaller as 
u — > oo. We find from (22) that the corresponding g(t) approaches the origin from the 
c*-direction. In order to show the existence of such solutions, we decompose M into a 
block diagonal form. Namely, 

™=( £ o F £ °oH A - < 23 > 

where eP is a k x k matrix whose eigenvalues have a negative real part. Similarly, eQ is 
an (n — k) x (n — k) matrix, where its eigenvalues have a positive real part. Note that 
P, Q and A are independent of e because M is linear in e. Define the tilde operation 

x = FT l x, X = R- 1 o X o R (24) 

for a point x G E and a map X : E —> R n , e.g., F(c) = R- l F(R~c) = R~ l F(c). The 
RGE (18) can be written as 

= eAb{u) + H{u, b(u)). (25) 

du 

If To prove (20), we use (5) again. 

V(c* +h) = V(c*) + DV(c*)h + o(|A»|), 

while the left-hand side is 

V(c* +h) = e 2 V(a* + h/e) = e 2 V{a*) + eDV (a*) h + o(\h\). 
Comparing the linear term in h, we have (20). 
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Let 



Then 



and 



ATT A r T 1 

4- = eAU(u), = eAT(u) (27) 

du du 



e eAu = U(u)+T(u). (28) 

We focus on a solution that behaves as b(u) — > as u — > oo. The integral equation 
corresponding to it is 

5( u ) = U(u)p + / du'l/(u - u')H(u', b{u')) - / du'T(u - u')H{u', b{u')), (29) 

JO Ju 

where p — (pi, ...,Pk, 0, 0) specifies an initial condition in the following way: 
&(0)i=pj for i = 

b(0)i = - (^j™ du'T(-u')H(u',~b(u'))^ for i — k + 1, ...,n. (30) 

We can show that (29) has a unique solution if e and p are sufficiently small. 
Moreover, we find that the solution satisfies 

b(u)\ < Je~ aeu (31) 

for some J > 0. Here a is a positive number, such that —a is strictly greater than 
the real part of every eigenvalue of P. In order to prove (31), we need to extend the 
usual proof of the stable manifold theorem (see [8]), which is applied to the case of an 
autonomous system, to the case of the nonautonomous system (18). Details of a proof 
of (31) are found in Appendix A. 



4. Universal Asymptotic Form of b(u) 

In this section, we derive a universal asymptotic form of b(u) by applying (31) to 
the right-hand side in (29). (Here, "universal" means that the asymptotic form is 
independent of p. ) 

For this purpose, we give a more concrete form of P. Let A; (/ = 1, ...,n_) be the 
distinct eigenvalues of P with the multiplicity d\. We denote by W\ the generalized 
eigenspace associated with A;. Clearly, dimW 7 ; = d\. Taking an appropriate basis for 
R n , P is represented as a block diagonal form. Here the I th block P t is a di x di upper 
triangle matrix whose diagonal components take a common value A/. Furthermore, the 
basis allows us to assume that iVj = Pi — Xild t is a nilpotent matrix, namely N Vl ~ x ^ 0, 
N Ul = for some 1 < v\ < d\. An arbitrary element x G R n can be decomposed as 

x = J2* {l) + x^eW h (32) 

i=i 
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where ajW is an element of the subspace spanned by the positive eigenmodes of M. 
Applying U (u) to the both sides, we have 

U{u) X = Y.^ lU Y. [j %T^* {l \ (33) 
1=1 k=o K - 

with the convention A z ° = 1^ even if Ni is the zero matrix. 

We can show that 4 " 

Mc* = -ec*, (34) 

which implies Pc* = — c*. Then we set 

Ai = -1. (35) 

Using (33) and (35), we can estimate the right-hand side of (29). The first term is 
written as 

U(u)p = Y,Gf\u)e^ u , (36) 
i=i 

where Gi\u) is a polynomial of degree at most v\ — 1. Since (36) explicitly depends on 
p, it is non-universal. 

In order to obtain a universal asymptotic form, universal terms should dominate 
over (36) when u — > oo. Let us find a condition that such terms appear from the 
remaining part. The integral containing U in (29) is divided as 

f'U f'U 

U(u-u')e- eu 'F(c*)du' + J o U{u-u') (h{u\ b{u')) - e~ eu ' F (c*)) du' . (37) 

Applying (33) to F(c*), the first integral is easily calculated. It is important to notice 
that the case of / = 1 has to be treated separately, because the factor exp(— eXiu') in 
U(u — u') cancels exp(— eu') in the integrand. When I = 1 and k — v\ — 1 in (33), 
the cancellation brings a term proportional to u Ul exp (—eu), which is not contained in 
Gf\u). Writing this explicitly, the integral is expressed as 

f U(u - u')e- eu 'F (c*) du' = (m) ^ - N?' 1 ^ (c*) + V e x ' eu G?(u). (38) 
Jo ev\\ 

Here, it is straightforward to check that G<p (u) is a polynomial whose degree is at most 
vi — 1, the same order as (36), so that its universal behavior is obscured by the non- 
universal nature of U(u)p. On the other hand, the first term can dominate over (36) in 
the case when 

3RA Z < — 1 (Z = 2,...,n_). (39) 

+ In order to prove (34), note that 

V (c* + he*) = V (c*) + DV (c*) he* + o(\h\). 
On the other hand, the left-hand side is equal to 

(l + h) 2 V(c*) = -e{l + h) 2 c* 
because of (5) and (12). Comparing the linear term in h, we get (34). 
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Moreover, as we see in Appendix B, this term is most dominant in the right-hand side 
of (29) when (39) holds. Thus, we obtain 



v i n —eu 



b(u) 



(eu) Ul e 
ev\\ 



N f-ip>W (5*) + O (u^- 1 ^) 



(40) 



Let us revert to the original variables by (22). We get, under the condition (39), 



m = 



et + 1 



ea 



: , e 3 (In (et + 1))- ATVi _ ipi i } \ + Q f (ln(rf + l))' 



e(et+l)u 1 \ 



(et + 1)' 



d* , (In , 



i/i-i' 



-N^F' ' (a*) + O 



t 2 



(41) 



as t — > oo. This is the main result of this paper. It is worthwhile to note that the (hit)" 1 
term appears, which brings about a (lnlnL)^ 1 correction in general. 



The case of dim W = 1 The simplest case is that dimlT^i = 1. In this case, using the 
unit eigenvector e* = a* / \a*\ in Wi, 

Ni =0 

vi = 1. (42) 

Then the result (41) is simplified to 



g (t) = f + l f(F(a*),e*)Re* + 0(±). 



(43) 



Note that we have removed ' ' from g and a* by applying R. Now let us write R as a 
set of column vectors 



R = v n ) . 

Similarly, we write in terms of a set of row vectors 

/ «i \ 

zr' = " 2 



(44) 



(45) 



\ u n J 

Since R~ l R — I n , we have 

(«i,t*j) = <Jij. (46) 

We know that a* is an eigenvector with the eigenvalue —1, so we can take t>! = a*. In 
this case, 

/ 1 \ 




a = R a 



V o y 



e*. 



(47) 
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Therefore 

=F 1 (a*), (48) 

where the right-hand side is the first component of F(a*). Using this convention, the 
asymptotic behavior (43) is simplified further to 

9(t) = ^ + ^F 1 (a*)a* + 0(±). (49) 





t t 2 1V ' \t 2 
5. Application to the two-dimensional XY model 

In this section, we illustrate our method using the two-dimensional classical XY model 
[9]. The beta function up to subleading order of this model and the two-point correlation 
function containing loglog correction are originally derived by Amit et al . [6]. They 
obtained the asymptotic form of the coupling constants by explicitly integrating the 
RGE. Here we rederive the asymptotic form within our formulation. 
The 2D classical XY model has the following RGE [6, 10]: 

d 9i 2 D 2 

-^T = ~92 ~ B ig 2 gi 

(50) 

, > - - , ~ B f? ) ■ (51) 
-9i92 J V ~ A i9 3 2 J 

Solving (14), we get a non-trivial solution 

"*=(!)■ ( 52 ) 
Inserting this into (21), one finds 

M =(-i ~o 2 )- < 53 > 

The eigenvalues and corresponding eigenvectors of M are 

-1 <-> Q and 2^ ^~ 2 ^. (54) 

Namely, the space of negative eigenmodes of M is one-dimensional. It indicates that 
the critical surface along a* is in fact a line. The transformation matrix R and 
the diagonalized matrix A are obtained from the eigenvectors and the eigenvalues 
respectively. The result is 
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The condition (39) is satisfied and dim W\ = 1 in this example. Furthermore, since we 
chose a* as the first column in R, (49) is applicable. Using R, we compute 

*(*•) =*-*(«•) = -i ( ^ ) . (56) 
From the first component of the above result, we conclude that 



a* 1 log t * _ / 1 . 

a* + O - (57) 



i 2 t 2 Vt : 
for the critical line. This is consistent with the original result 



6. Summary and Discussion 

We have obtained an asymptotic formula for multiple marginally irrelevant coupling 
constants in the case where the two-loop beta functions are known. We first change the 
variables in the given RGE as (6) and (8). One can extract a linear part M defined in 
(19) from the resulting differential equation if there is a real solution c* for (12), although 
the original RGE cannot be linearized. It is assumed that there are no eigenvalues of 
M with zero real part in the present investigation. 

Next we have shown that, if we take e in (8) to be sufficiently small, there is a 
k dimensional neighborhood iV of c* such that trajectories of the RGE starting in N 
approach the origin along c*, where k is the number of eigenvalues of M with a negative 
real part. Furthermore, if the eigenvalues of M with a negative real part satisfy the 
condition (39), the asymptotic formula of g(t) becomes universal and is given by (41). 

As we have concretely shown by using the two-dimensional classical XY model, the 
advantage of this formula is that we do not need to integrate the RGE explicitly. 

The non-linearity having the original RGE is changed into solving (12) in our 
formalism. All quantities appearing in (41) can be computed using simple linear algebra. 
Thus, our approach can be applicable even though the RGE is too complicated to 
integrate. Application to such a complicated RGE is a future problem. 

Finally, we comment on our previous investigation on this topic. As we mentioned 
briefly in the introduction, another RGE is derived with respect to the scaling invariance 
possessed by the leading-order RGE in reference [2], where the asymptotic formula in 
the in the lowest order is obtained. It is consistent with the present work. The rela- 
tionship between the linear part of the new RGE and M in (19) is clarified in reference 
[3]. However, a direct relationship between M and the original RGE was not clear. In 
this paper, it is found that M naturally comes out in the RGE by the change of vari- 
ables, which makes it possible to derive the asymptotic formula for a beta function up 
to the subleading order. The previous formulation can also be applied to deriving the 
correlation-length exponent of phase transitions in infinite order, which will be extended 
to the case of the higher-order beta function. 
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Appendix A. Existence of the solution for (29) 

The purpose of this appendix is to show existence of a unique solution for the integral 
equation (29) and derive the estimation (31). To this end, we need estimations of H, 
U and T. Throughout this appendix, we omit "~" for brevity, (e.g., we write b or H 
instead of b or H respectively. ) 

Appendix A. 1. Lipschitz-type condition for v and f 

We first derive a Lipschitz-type condition for v and / defined through the following 
equation: 

V (c* + b) =V (c*) + DV (c*) b + v(b) (A.l) 
F(c* + b) = F(a*) + DF(a*)b + f(b), (A.2) 

which is employed for the estimation of H. Let 

N v = {b;beR n , \b\ < v }. (A.3) 

First we prove that: for any £1,^2 > 0, there exists a number 77 > such that 

|«(bi)-«(b2)| <Zie\b 1 -b 2 \ (A.4) 
1/(60-/(62)1 <&\b 1 -b 2 \ (A.5) 

for &!, 6 2 £ N er] and < e < 1. 

Proof. We define v by the following equation: 

V (a* + b') = V (a*) + DV (a*) b' + w (b')- (A.6) 
For b' 1; b' 2 , let s = b' 2 — b' x and define 

= (b; + 0a) . (A.7) 

Note that V(0) = Vo(b' 1 ) and Y(l) = v (b' 2 ). Taking the derivative with respect to 9, 
we have 

= Dv (b[ + 9s) s. (A.8) 
Integrating from to 1, we get 

\v (b[) - v (b' 2 )| = |Y(1) - F(0)| < f 1 d9 \Dv (b[ + 9s)\ \s\ . (A.9) 

Jo 

Since Dv is continuous and Dv (0) = 0, for every ^1 > there exists a number rj > 
such that x e N v implies \Dv (x)\ < £1. For bi,b 2 G A^, since bi + 6>s G N v , (A.9) 
leads to 

|« (6i)-»o(^)|<6|6i-^l- (A.10) 
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Next we consider estimation in a neighborhood of c*. Putting b = eb' and using 
the scaling property of V, we have 

V (c* + b) = e 2 V (a* + b') 

= e 2 (V (a*) + DV (a*) b' + v (b')) (A.ll) 
Comparing (A.l) and (A.ll), we have 

v(eb')=e 2 v (b'), (A.12) 
which leads to 

\v(b!) - v(b 2 )\ = e 2 \v (b[)-v (b 2 )\ 



<e 2 Z 1 \b' 1 -b' 2 \ 

= e£ 1 |b 1 -6 2 |, (A.13) 
for 61, b 2 £ A err Repeating a similar argument for /, we get 

- f(b 2 )\ <e%\b 1 -b 2 \, (A.14) 
which implies (A. 5) because < e < 1. 

Appendix A. 2. Lipschitz-type condition for H(u, b) 

Using (A. 4) and (A. 5), we readily obtain a Lipschitz-type condition for H in the 
following form: for any £ > 0, there exists i] > such that 

\H(u, b ± ) - H(u, b 2 )\ < (£e + we 2 ) |&i - b 2 \ (A.15) 

for all 61, 6 2 G A er? and < e < 1. Here, 

w = \DF(a*)\. (A.16) 

Proo/ Recall 

H(u,b) = v(b) + e- €U F{c* + b). (A.17) 

Then 

\H(u,b 1 )-H(u,b 2 )\ < | V (6 1 )-«(6 2 )|+e- eu | J DF(c*)(b 1 -62)| 

+ 6—1/(60-/(62)1 
< \v (bi) - v (b 2 )\ + e 2 w \b ± - b 2 \ 

+ |/(6i)-/ (6 2 )|, (A.18) 

where we have used 

|DF(c*)| = e 2 u>. (A.19) 

Using (A.4) and (A.5) for 6 = 6 = £A we get (A.15). 
Corollary: if (A.15) holds, then 

\H (u, 5i)| < (6 + e 2 u?) |bx| + e - ae V |F (a*) I , (A.20) 

where < a < 1. 

Proof. Setting 6 2 = in (A.15), we immediately get 

H (u, 61) - e~ eu F (c*)| < (6 + e 2 w) |6i| , (A.21) 
which implies (A.20). 
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Appendix A. 3. Estimation of U andT 

Following reference [8] (Section 2.7), we derive: there are positive constants a, a and 
K, such that 

\T(u)\ < Ke utu (forw<0) 

\U(u)\ < Ke- {a+a)eu (ior u > 0), (A.22) 
where < a < 1. 

Proof. Let fij (j = 1, ■■■,n + ) be distinct eigenvalues of Q. Here 9?/ij > for all j. We 
can choose sufficiently small o and sufficiently large K such that 

5fytj > a > for all j = 1, n + 

\T{u)\ < Ke utu for all u < 0. (A.23) 

Similarly, let A; (/ = l,...,n_) be distinct eigenvalues of P. We take positive a' 
satisfying 

3^ < -a' < (A.24) 

for all / = 1, ...,n_. Then there exists K' > such that 

|C/(u)| < ir'e" '™ (A.25) 

for all u > 0. It should be noted that a' < 1 because of (35). Further, if we choose a so 
small that 

9^ < -a' < -a < 0, (A.26) 
then we can write 

a' = a + a (A.27) 
using some < a < 1. Combining (A.23) and (A.25), we obtain (A.22). 

Appendix A. 4- Existence of a solution 

We define the space C of continuous mappings on [0, oo] into R n in the following way: 
C = [b; \b(u)\ < er]e~ aeu , u > o} . (A.28) 

Define the metric on C by 

p (&!, b 2 ) = sup (u) - b 2 (it)) e aeu \ , (A.29) 

for bi, b 2 e C. Note that if b e C, b(u) G N erj for all u > 0. Then according to (A. 15), 
for £ > 0, there exists rj > such that 

< (te + we 2 )p(b 1 ,b 2 )e- aeu (A.30) 

holds for bi, b 2 e C, w > 0. 
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Next we introduce the following mapping ip on C: 

ru roo 

$ (b) (u) = U{u)p + / du'U(u - u')H(u', b(u')) - / du'T{u - u')H(u', b(u')).(AM) 

JO Ju 

We have 

|V (61) («) - ^ (62) («)| < # (£e + ew 2 ) p (bx, b 2 ) 

x ^ e - (a+<T)e{M - M,) e- w dM' 

+ e^ u ~ u '^e- aeu 'du'^ , (A.32) 
according to (A. 22) and (A. 30). Here the first integration is evaluated as 

ru 1 ~—aeu 

/ e -(a-Hr) e (u-u') e -a eU ' dM / = _ (g-rcu _ e -(a+a)e U \ < _ /^ggN 

Jo ere v ' ere ' 

The second one is 

/ e CT(u - u ' )e e- aeu 'dw' = -e~ aeu < . (A.34) 

Ju (a + er)e ere 

In this way, we obtain 

2K 

p (60 , ij (6 2 )) < — (£ + cu;) p (61, 62) , (A.35) 
cr 

for b 1 ,b 2 e C. 

Now we show that ip is a contraction map on C into itself if we choose £, 77, e and 
p such that the following inequalities hold: 

K (\p\ + \F (a*)|) <^ (A.36) 

— (£ + e^) <1 (A.37) 

cr z 

It is easy to see that such choice is in fact possible.* 

When (A.36) and (A.37) hold, we can show that ip maps C into itself. Let us recall 

(A.20). If b e C, then 

\H (u, b(u))\ < (en (f e + e 2 w) + e 3 \F (a*)\) e~ aeu . (A.38) 

* For example, we first fix £ = £0 such that < £0 < °7(4-?0, which is realized by an appropriate 
choice of r)(— rj ) according to (A. 15). Next we choose e = e satisfying 



eo < mm 



(iK \F (a*)\' w (iK * )) " 



By construction, (eo,£o) satisfies (A.37). Furthermore, choosing p such that 

° <|p|< ; W\F(a*)\ - 

we find that (A.36) holds. 
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Combining this and the estimations (A. 33) and (A. 34), (A. 31) yields 

\ip (6) (u)\ < Ke~^ a+ ^ eu \p\ + (eTy (f + ew) + e 2 \F (o*)|) — e""™ 

/ IK \ "IK 

<K(\p\ + —e 2 \F (a*) | J e" aeu + er] (f + ew) — e~ aeu 

< er]e- aeu (A.39) 

for all u > 0. In the last inequality, we have used (A. 36) and (A. 37). It shows that 
ip (b) E C. Furthermore, (A. 35) and (A. 37) indicate that ip is a contraction map. 
Thus, existence b e C of the solution 

b = i/j(b) (A.40) 

and its uniqueness follow from the fixed point theorem (see the textbook [11], for 
example). Since b belongs to C, we find that there is a positive number J such that 

\b(u)\ < Je~ aeu (A.41) 

for u > 0. 

Appendix B. Estimation of the remaining terms in (29) 

According to (29) and (36)-(38), it is sufficient to show that 

ru n - 

/ U(u — u') (H(u', b(u')) — e~ tu F (c*)) du < ^ G?(u)e Ai£U (B.l) 



/ T{u-u')H{u',b{u'))du' 

Ju 



< Be- 



rn) 



for all u > 0, in order to complete the derivation of (40). Here B is a positive constant 
and G\p(u) is a polynomial whose degree is at most v% — 1. 

Let us first show (B.l). According to (A. 2), H is written as 

H (u, b(u)) = v (b(u)) + e~ eu (F (c*) + (c*) + / (&(«))) . (B.3) 

Since V, F e C 2 (E), (j the reminders v, / also belong to C 2 (E). Using the Taylor 
theorem, we find that there exists a number (0 < 9 < 1) such that 



v ' jk 2dbjdb k 



(B.4) 



holds. Employing the estimation (31) for b(u), we can show that for sufficiently large 

-2aeu 



v (b(uj) \ < B lC - 
A similar estimation holds true for /, i.e., for some B 2 > 0, 
~f (b(u))\ < B 2 e 



-2aeu 



(B.5) 
(B.6) 



sec the second footnote in Section 2. 
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If the assumption (39) holds, we can take 



- < a < 1. 
2 



(B.7) 



Further, we choose a such that 2a ^ —A; for all / = 1, ...,n_ for later use. Combining 
the condition (B.7) with (B.5) and (B.6), it follows that there is a number B 3 > such 
that 



H(u,b(u)) -e~ tu F{t*)\ < \v (b(u))\ + e~ eu ( 



DF (c*) b(u) + f(b(u] 

< B x e~ 2aeu + DF (c*) e - (1+Q)eu + B2e ( - 1_2a)< 

< B 3 e~ 2aeu 



(B.8) 



for all u > 0. Since the same estimation holds for every component projected onto W h 
we obtain 



\u(u - u') (H(u', b(u')) - e~ £U F (c* 



1=1 



k=0 



k\ 



H K ''(u', b(u'))-e- eu F {l) (c*) 



i/j-i 



J=l fc=0 ^' 

Integrating the both sides by u' over [0,-u], we obtain (B.l) 
As for (B.2), we use (B.8). Then we have 



(B.9) 



H (V, b(u')) | < |.F(c~*)| e- eu ' + B 3 e~ 2aeu ' < Ce 



(B.10) 



for some C > because we have set a > 1/2. Then the integral is easily estimated as 



du'T{u - u')H (V, b(u' 
Thus we get (B.2). 



< KC / du'e ae{u - u,) e- eu ' 

J u 



KC 



(<r + l)e 



(B.H) 
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